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Abstract 

We calculated the matter effect, including both the Earth and supernova, on the detection of 
neutrinos from type II supernovae at the proposed Daya Bay reactor neutrino experiment. It is 
found that apart from the dependence on the flip probability Pjj inside the supernova and the 
mass hierarchy of neutrinos, the amount of the Earth matter effect depends on the direction of the 
, incoming supernova neutrinos, and reaches the biggest value when the incident angle of neutrinos 

is around 93°. In the reaction channel v e + p — ► e + + n the Earth matter effect can be as big as 
about 12%. For other detection processes the amount of the Earth matter effect is a few per cent. 
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I. INTRODUCTION 



One of the most spectacular cosmic events is a supernova (SN) explosion which ends the 
existence of a giant star. There are two types of such explosions: types I and II. In the case 
of type I supernovae the explosion gives out gigantic firework display, while in the type II 
supernovae intensive neutrinos are emitted for a very short period of time followed by an 
intensive electromagnetic radiation. Supernovae are a very good natural laboratory to study 
various fundamental issues of physics and astrophysics. 

The type II SN explosion is caused by the mechanism of core collapse. The total energy 
release is approximately the gravitational binding energy of the core. It generates intensive 
neutrinos which take away about 99% of the total energy. The explosion itself consumes 
about 1% of the total energy. The neutrino emission happens several hours prior to the 
brightening of the SN which only consumes 0.1% of the total energy. The observation of the 
SN neutrinos can serve as an early warning for the optical emission of a type II SN. 

The first observations of the type II SN were made by Kamiokande II , 1MB, and Baksan 
in 1987, with totally 24 neutrino events observed in a short period of about 13 seconds 
PJ. Although this number of events is too small for a quantitative study of neutrinos from 
supernova explosion, they are definitely very valuable for providing the first signal of cosmic 
neutrinos and initiating astrophysics and particle physics studies of supernova neutrinos. 

The Daya Bay underground neutrino laboratory has been proposed and R&D is now 
being carried out. It is located in south China close to Hong Kong and provides an excellent 
site for measuring the neutrino mixing angle #13 by using 4 reactor cores in two clusters, each 
core with a thermal power of 2.9 GW and therefore 11.6 GW total power. Furthermore, two 
more cores with an additional thermal power of 5.8 GW are expected to be online in 2011. 
Since the site is close to a group of hills, the cosmic ray background to the neutrino signals 
is vastly reduced. While the main purpose of the Daya Bay experiment is to measure a 
neutrino mixing angle, it can be used to detect possible neutrino events from a SN explosion 
and serves as a part of the Supernovae Early Warning System (SNEWS) 0. Hence the 
theoreticalprediction for the detection of SN neutrinos at Daya Bay is very desirable. 

In Ref. |3J] the authors calculated the expected neutrino events from a type II SN explosion 
in Borexino. For a typical SN at a distance of 10 kpc a burst of around 110 events would 
appear in the detector of Borexino. In this calculation, matter effects are ignored. Later 
Whisnant and Young 0] considered the detection of type II SN neutrinos in the Daya Bay 
experiment with the matter effect inside the supernova being included. This leads to the 
modified results which depend on the flip probability of neutrinos while passing through 
resonance regions inside a supernova. 

In the realistic case, the neutrinos produced from a SN explosion will likely go through 
some portion of the Earth before reaching the detector. Therefore, in order to give a more 
accurate prediction for the number of signal events, the Earth matter effect should also be 
taken into account. The purpose of the present paper is to calculate the Earth matter effect 
and therefore takes into account all matter effects, both the Earth and Supernova, on the 
detection of type II SN neutrinos in the Daya Bay experiment. 

II. THE FORMALISM 

During a SN explosion, neutrinos are produced in two bursts. In the first burst which 
lasts for only a few milliseconds, electron neutrinos are generated via the inverse /5-decay 



process e~ + p — > v e + n which leads to a neutron rich star. In the second burst which is 
longer (O(10) seconds) neutrinos of all flavors (u a and v a with a being e, /i, r) are produced 
via the e + e~ annihilation e + + e~ — > i/ a + i> a , the electron neutrino antineutrino annihilation 
u e + z/ e — > z/ Q + z/ Q , and nucleon-nucleon bremsstrahlung N + N — > N + N + v a + v a processes. 
While these neutrinos generated in the SN explosion travel out of the SN, they interact with 
the high density core of the SN. Neutrinos v x (v x = v^v T ^v^v T ) decouple first from deep 
inside the core due to their lesser interactions as they only interact via the neutral current. 
Then v e decouples. The last one to decouple is v e due to its interaction with the neutron 
which has a higher density than the proton. The typical temperatures of these neutrinos 
are in the following ranges 0: 

T, e = 3-4MeV, T Pe = 5 - 6 MeV, T Vx = 7 - 9 MeV. (1) 

For the neutrino of flavor a, the time-integrated neutrino energy spectrum can be de- 
scribed by the Fermi-Dirac distribution, 

7V(°) E 2 

F<-°)(E) = — - (2) 

a 1 1 F 2 T*exp(E/T a ) + V 1 } 

where E is the energy of the neutrino and T a the temperature as given in Eq. ([T|) . iV^ ) is 
the total number of the neutrino of flavor a, and Fj, where j is an integer, is defined by 



Fj = / dx. (3) 

Jo exp(x) + 1 

The average neutrino energy is then obtained from the distribution function Eq. (j2J) 

<^ 0) > = W 

The total number of the neutrino u a is 

where is the luminosity of u a . Then the energy spectrum function can be rewritten as 

E 2 

fM(E) = . (6) 

Q { ' F 3 T* exp(E/T a ) + 1 1 ) 

In the simplest argument the luminosity of each flavor is the same 0,0]. Hence 

r(0) _ °-"p(0) m 

where Egpf is the total energy released during the SN explosion. 

When the SN neutrinos of each flavor are produced they are approximately the effective 
mass eigenstates due to the extremely high matter density environment. While they propa- 
gate outward to the surface of the SN they could experience level crossing (neutrinos jump 
from one mass eigenstate to another) in the MSW resonance regions 0. There are two MSW 



resonance regions which are determined by the two pairs of parameters (Am^, sin 2 20i3) 
which is referred to as the high resonance region and (Amjj, sin 2 26 l i2) the low resonance 
region, where Am^ = m\ — m 2 , i.e., the mass square difference of the k and j neutrinos, 
and #jk the corresponding mixing angle. Let us denote the probability that the neutrinos 
jump from one mass eigenstate to another at the high (low) resonance layer by Ph{Pl)- The 
large mixing angle (LMA) solution of the solar neutrino constrains the flip probability at 
the low resonance region to be zero (Pl = 0) 0. Ph is a function of the mixing angle sin#i3 
and Am| x . Since the value of the mixing angle sin 2 2#i3 is still not determined, in numerical 
calculation we let Ph vary between and 1. The level crossing diagrams are different for 
normal and inverted mass hierarchies [5|. This leads to different forms for neutrino flux of 
mass engenstates at the surface of SN. 

When the neutrinos arrive at the Earth after travelling through the cosmic distance all 
the oscillation factors are averaged out and there is no coherence between different mass 
eigenstates. Hence the neutrino fluxes reaching the Earth are the following: For the normal 
hierarchy (Am^ > 0), 

Fg° = f£ + (\U*\ 2 Ph + \UM1-Ph)){F$!>-FW), 

Ft = Fff + \UM^-F^), 
2Fjf ) = F® + F® - (\U e2 \*P H + |t/ e3 | 2 (l - P H )){F^ - F(% 

2Ft = 4 0) + F® - \U el \\F^ ~ F®), (8) 

and for the inverted hierarchy (Am| x < 0), 

F%> = F® + \U0\*(FW-FW), 

F£ = FW + (\U el \ 2 P H + \U e3 \\l - P B ))(F® - F(% 
2F (i) = F (o) + F (o) _ |f/e2|2(F (0) _ F (o })) 

2F« = pf + F«?> - (\U el \*P H + |^ e3 | 2 (l - Ph))(F^ - Ff), (9) 

where U e i(i = 1, 2, 3) are the elements of the neutrino mixing matrix. 

Due to the smallness of the mixing angle #13, we ignore the terms proportional to sin# 13 , 
then Eqs. (jSJ) and can be simplified as follows: For the normal hierarchy, 

Fjf) = P H ^6 l2 F^ + (1 - P H sm% 2 )F^\ 

pjf } = cos% 2 F^+sm% 2 F^\ 
2FW = (l-P H sm% 2 )F^ + (l + P H sm 2 6 12 )F£\ 
2Fjf ] = sm% 2 F^ + (1 + cos 2 ^ 12 )Fi°), (10) 
and for the inverted hierarchy (Am^ < 0), 

F« = sin 2 ^ 12 F( e °)+cos 2 ^ 12 FW, 

FiP = P H cos% 2 F^ + (1-P hC os 2 9 12 )fI / J, 

2F« = cos 2 fl 12 Pf + (1 + sin 2 £ 12 )Fj°), 

2Fg> = (1 - P„cos 2 # 12 )Pj e 0) + (1 + Phcos 2 6 12 )f£\ (11) 

If the SN matter effect is neglected, the neutrinos, which are produced in the core of 
the SN and propagate through cosmic distance to reach Earth, are subject to the vacuum 



oscillation only. In this case when the Earth matter effect is also neglected the neutrino 
fluxes at the detector are given by: 

FSP = (l-isin 2 2^ 12 )F( e ) + isin 2 2^ 12 FW, 

F$? = (l-^sin 2 2^ 12 )Ff + isin 2 2^ 12 F(°), 

2F$P = ^sin 2 2^ 12 FW + (2-isin 2 2^ 2 )FW, 

2FZ ] = \^29 l2 F^ + (2 - \^29 l2 )F^. (12) 




FIG. 1: Illustration of the path of the SN neutrinos reaching the detector in the Earth. D is the 
location of the detector, 9 is the incident angle of neutrinos, and O is the centre of the Earth. 



In Eqs. (jlOjl and ([TTf the Earth matter effect is neglected. In reality, neutrinos from an 
SN may likely travel through a significant portion of the Earth before reaching the detector 
and are therefore subject to the matter effect. The Earth matter effect can be parameterized 
in terms of the distance between the neutrino and the Earth center and the total distance 
travelled in the Earth before reaching the detector. The distance that the neutrino travels 
in the Earth depends on their incident orientation. Suppose a neutrino reaches the detector 
with the incident angle 9 (see Fig. 1), then the distance the neutrino travels through the 
Earth is 

L = (-R + h)cos9 + ^R 2 -(R- h) 2 sm 2 9, (13) 

where h is the depth of the detector in and R the radius of the Earth. The distance of 
the neutrino to the Earth center can be obtained as follows: Let x be the distance that the 
neutrino travels into the Earth, which is shown as the segment AB in Fig. 1. The distance 
of the neutrino to the center of the Earth x is given by 



x = sJ(R - h) 2 + (L- x) 2 + 2(R -h)(L- x)cos9. 



(14) 



The Earth matter effects on the detection of SN neutrinos have been discussed extensively 
in literature (H, EU EH EH, EH Eil E3 As mentioned before, after travelling the cosmic 
distance to reach the Earth, the arriving neutrinos are definite mass eigenstates which, 
then, oscillate in going through the Earth matter. Let Pj e be the probability that a neutrino 
mass eigenstate Vi enters the surface of the Earth and arrives at the detector as an electron 
neutrino u e . Then the flux of u e at the detector, denoted as F®, can be written as 

*2 = £flA (is) 

i 

where Fi is the flux of V{ at the Earth surface, in either the normal or inverted hierarchy. Pj e 
obey the unitary condition Ya Pie = 1. In the case where the Earth matter effect is ignored, 
Pie = \U e i\ 2 , then one recovers the expressions of Eqs. ([%])- (fTTj) . After some straightforward 
derivations with the terms proportional to sin# 13 being ignored, one obtains, in the case of 
the normal hierarchy, 

f d W = p 2e P HF (o) + (1 _ P 2e P H ) F £\ 
F°} N) = (l-P 2e )F^ + P 2e F^\ 
2F°W = (l-P 2e P H )FM + (l + P 2e P H )F£\ 

2F°} N) = P 2e pf + (2 - P 2e )pj^ (16) 
and in the case of the inverted hierarchy, 

= P 2e Pi 0) + (l-P 2e )Pi 0) , 

f£ {i) = p h (i - p 2e )Pi e 0) + (l + p 2£ p h - p h )f£\ 

2F^ = (l-P 2e )Pi e °) + (l + P 2e )pW, 

2P£ (/) = (1 + P 2e P H - P H )F$? + (1+P H - P 2 ePu)F^. (17) 

The probability p e (i = 1,2,3) has been calculated in Refs. Q and [l(j. In the Earth 
the potential which an electron neutrino experiences is 

V(x) = V2G F N e (x), (18) 

where Gf is the Fermi constant and N e (x) is the Earth electron number density. Let us 
introduce a small parameter 

, . 2EV(x) 
ex = \' . 19 

For a typical neutrino energy E = 10 MeV, e _is less than 0.13. Neglecting contributions of 
0(e 2 ), P 2e can be expressed as the following 



P 2e = sin 2 fl 12 + J-sin 2 2# 12 H dxV (x)sm<j)™ (20) 
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where Q\ 2 ~ 33° [16j and 4>™_^ b is defined as 



dxA m (x), (21) 



where 

A rri ^ i 

A m {x) = ^(cos 29 12 - t{x)Y + sin 2 20 12 . (22) 

Let the matter density inside the Earth be p(x) then the nucleon number density is 
2p(x)/(m p + m n ), which, for matters of equal number of protons and neutrons, leads to the 
following density of electrons inside the Earth: 

N e (x) = p(x)/(m p + m n ). (23) 

In our calculations of the Earth matter effect we choose a simple model for p(x) which 
is called the mantle-core- mantle density profile given in |l7[ and |18j . In this model a 
step function is used to describe the Earth matter effect: p = Ylgjcrn? for the core, and 
p = 5g/cm 3 for the mantle. The core radius and the thickness of the mantle are each half 
of the Earth radius. 

In the next section we will use the formulas given in this section to calculate numerically 
the neutrino fluxes at the detector, taking into account the Earth and supernova matter 
effects, and then apply the results to various physical channels through which neutrinos can 
be detected. 



III. NUMERICAL RESULTS 

In the following we calculate the predicted numbers of events that can be observed through 
various reaction channels at the Daya Bay experiment. This will be done by integrating, 
over the neutrino energy E v , the product of the target number Nt, the cross section of each 
channel, and the neutrino flux function F^(E u )/AttD 2 given in Eqs. (Jini) and dnj) @: 

N(i) = N T J dE v a{i){Ev)^ 2 F°{E v ), (24) 

where D is the distance between the SN and the Earth and the index i represents different 
channels through which the SN neutrinos are observed. 

The target material includes electrons, protons and carbon nuclei since the chemical 
composition of the expected liquid scintillator, mesitylene or pseudocumene, is C 9 H 12 which 
is the same as that of Borexino 0. The cross section for each reaction channel can be 
found in Ref. For neutrino-electron elastic scatterings via charged and neutral currents, 
v a {pu) + e~ — > Voipa) + e~, the cross sections are: 

<r(u e e) = 9.2 x 10" 45 ^(MeV) cm 2 , 

a{p e e) = 3.83 x 10" 45 ^(MeV) cm 2 , 

a{y^ T e) = 1.57 x 10" 45 ^(MeV) cm 2 , 

a(u^ T e) = 1.29 x 10~ 45 ^(MeV) cm 2 . (25) 

The inverse beta decay u e + p — » e + + n has the most number of events due to its large 
cross section and low threshold of the neutrino energy. 



a(P e p) = 9.5 x 10- 44 (^(MeV) - 1.29) 2 cm 2 , 
E th = 1.80 MeV. 



(26) 



For the neutrinos and 12 C system, there are two charged-current and three neutral-current 
reactions. Their effective cross sections are obtained by scaling the experimentally measured 
energy values from muon decay at rest to the energy scale for supernova neutrinos. In this 
way, one obtains the following effective cross section : 

charged-current capture of v e : v e + 12 C ^ 12 B + e + , 12 B ~^ 12 C + e~ + u e , 

(a( 12 C(v e , e+) 12 B)} = 1.87 x 10~ 42 cm 2 ; (27) 

charged-current capture of u e : u e + 12 C — > 12 N + e~, 12 N — > 12 C + e + + v e , 

{a( 12 C(u e , e-) 12 N)} = 1.85 x 10~ 43 cm 2 ; (28) 



neutral-current inelastic scattering of u, 

(a(u e l2 C)) 
(a(v e 12 C)> 



or v a : v + 12 C ^ 12 C* + v', l2 C* ^ 12 C + 7, 

= 1.33 x 1(T 43 cm 2 ; (29) 
= 6.88 x 1(T 43 cm 2 ; (30) 
= 3.73 x 1(T 42 cm 2 , x = //, r. (31) 



These reactions will be referred to collectively as neutrino-carbon scatterings. 

There will be several detectors located in the near and far sites. We take the total detector 
mass to be 300 tons which is the same as that of Borexino 0. Then the total numbers of 
target electrons, protons, and l2 C are 

= 9.94 x 10 31 , = 1.82 x 10 31 , nP = 1.36 x 10 31 . 

In the numerical calculations we take the total energy release of SN as 

E§1 = 3 x 10 53 erg = 1.97 x 10 59 MeV, 

and the distance of SN as 

D — 10 kpc = 3.09 x 10 22 cm. 
Other parameters are taken as follows: 

/i = 0.4 km, i2 = 6400 km, Am 2 21 = 7.1 x 10" 5 eV 2 , 9 12 = 32.5°. 

The numerical results for the event numbers for various reaction channels are obtained 
from Eqs. (|24| ) -(pT |) . (jTUj) and ()17|1. The results are shown in Figs. 2, 3, and 4. It can be 
seen from these figures that the Earth matter effect depends on the incident angle of the 
neutrino, the mass hierarchy, and the flip probability Ph- When the incident angle of the SN 
neutrino is smaller than a value 8q, where 6$ ~ 90°, the Earth matter effect can be ignored 
for all the relevant reactions that concern us here. 

The inverse beta decay u e + p — >• e + + n has the largest Earth matter effect among all the 
channels. The maximum Earth matter effect for the inverse beta decay appears at around 
6 ~ 93.6° in the cases of Pjj = 0, 1 (normal hierarchy) and Pjj = 1 (inverted hierarchy) 
and the effect is 12.0%. When the incident angle become larger than about 106°, the Earth 
matter effect becomes about 6.7% and almost independent of the incident angle. 



For the neutrino-electron and neutrino-carbon scattering channels the behavior is quite 
similar to that of the inverse beta decay. For neutrino-electron elastic scattering, the maxi- 
mum Earth matter effect appears at 9 ~ 92.7° and the amount could be as large as 3.75% 
in the case of Ph = 1 for both normal and inverted hierarchies, 3.11% in the case of Ph = 
for the inverted hierarchy, and 0.66% in the case of Ph = for the normal hierarchy. When 
the incident angle become larger than about 100°, the Earth matter effect is independent 
of the incident angle and the amount is about 2.16% for Ph = 1 (normal and inverted hier- 
archies), 1.76% for Ph = (inverted hierarchy), and 0.38% for Ph = (normal hierarchy), 
respectively. 

It should be noted that there are complications in dealing with reactions of neutrinos 
with 12 C . The effective cross sections in Eqs. (|27H3H) are given for supernova neutrinos 
without oscillations. It is pointed out in Ref. [3J| that the effective cross sections are affected 
when neutrino oscillations are taken into account. For instance, the oscillation of higher 
energy z/ M into v e results in an increased event rate since the expected v e energies are just 
at or below the charged-current reaction threshold. This leads to an increase by a factor of 
35 for the cross section l2 C(v e e~) 12 N if one average over a v e distribution with T = 8 MeV 
rather than 3.5 MeV. After numerical calculations for the reactions of neutrinos with 12 C, 
we find that the Earth matter effects become the largest when 9 ~ 91.8° and the maximum 
amount is 4.00% in the case of Ph = 1 for both normal and inverted hierarchies, 1.65% 
in the case of Ph = for the inverted hierarchy, and 2.00% in the case of Ph = for the 
normal hierarchy. Again when the incident angle become larger than about 100°, the Earth 
matter effect is almost independent of the incident angle and the amount is about 2.87% 
for Ph = 1 (normal and inverted hierarchies), 1.21% for Ph = (inverted hierarchy), and 
1.36% for Ph = (normal hierarchy), respectively. 

The above results can be understood from the fact that the oscillation behavior is deter- 
mined by the phase factor Am2 1 (eV 2 )L(m)/i?(MeV), where L is the distance in the Earth 
that the neutrino travels. Equation (|13|) gives the distance the SN neutrinos travel in the 
Earth. It can be seen from this equation that when 9 < 9 this distance is smaller than 
10 km and hence the amount of Earth matter effect is very small. In this case in Eq- 
(JUJ) and F? in Eq. (HHJ) can be replaced by their values in the vacuum \U e 2\ 2 and F Ve , 
respectively. When 9 is larger than 90° this distance becomes bigger than 100 km, the Earth 
matter effect becomes large and reaches a maximum for 9 ~ 91° — 94°. When 9 is larger 
than 100°, the distance that neutrinos travel in the Earth is greater than 2000 km, then 
there could be many oscillations in F®. This leads to averaging the Earth matter effect 
which is smaller than the maximum value. As an example, the flux of the SN neutrino v e 
as a function of neutrino energy in various cases are plotted in Fig. 5. It can be seen from 
Figs. 2 to 5 that the event numbers detected at the detector depend on the incident angle of 
neutrinos, the SN matter effects (Ph), and the mass hierarchy of neutrinos. By measuring 
these event numbers we can obtain information on the SN and the Earth matter effect, and 
the neutrino mass hierarchy. 




9 (in units of degrees) 

FIG. 2: The event number of the reaction i> e + p->e + + nasa function of the incident angle of the 
SN neutrino. The solid curve corresponds to Ph = 0, 1 (normal hierarchy) and Pjj = 1 (inverted 
hierarchy). The dashed curve corresponds to Ph = (inverted hierarchy). 



IV. SUMMARY 



We have calculated matter effects, including both the supernova and the Earth, on the 
detection of type II supernovae neutrinos at the proposed Daya Bay experiment. It is 
found that the amount of the matter effect depends on the neutrino incident angle, the 
neutrino mass hierarchy, and the flip probability Ph- When the incident angle of the SN 
neutrinos is smaller than a definite value # , the Earth matter effect is negligible due to the 
small distance that the neutrino traverses in the Earth. The Earth matter effect becomes 
manifested and reachs a maximum value for 9 ~ 91° — 94°. When 9 is greater than 100°, 
the Earth matter effect becomes insensitive to the neutrino incident angle and the effect is 
smaller than the maximum value. There is a window in 9, 9 ~ 90° — 100°, corresponding to 
the neutrino travelling distance of 500 km to 2000 km, in which the neutrino event numbers 
can vary significantly as a function of 9. This phenomenon can be understood by the fact 
that the Earth matter effect is controlled by Am^ and the variation of the neutrino number 
is significant for Aml 1 (eV 2 )L(m) / E v (MeV) ~ a few. 

In the reaction channel v e + p — > e + + n the amount of the Earth matter effect is the 
largest and the maximum amount is about 12%. For the neutrino-electron and neutrino- 
carbon scattering channels, the amount of the Earth matter effect is a few per cent at most. 
By measuring the event numbers in various channels we can obtain information on the Earth 
matter effect, the matter effect inside the SN, and the neutrino mass hierarchy. 
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FIG. 3: The total event number of neutrino-electron elastic scattering via charged and neutral 
currents u a (u a ) + e~ — > v a [p a ) + e~ as a function of the incident angle of the SN neutrinos. The 
solid curve corresponds to Pjj = 1 (normal and inverted hierarchies), the dashed curve corresponds 
to Pjj = (normal hierarchy), and the dotted curve corresponds to Ph = (inverted hierarchy), 
respectively. 
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FIG. 5: The flux of the SN v e at the detector multiplied by 10~ 8 . The solid curve represents the 
cases where supernova and Earth matter effects are included for 9 = 180°, Pjj = 0,1 (normal 
hierarchy) or Pjj = (inverted hierarchy). The dashed curve represents the case where supernova 
and Earth matter effects are included for 9 = 180°, Pjj = (inverted hierarchy) or the case where 
the supernova matter effect is included but the Earth matter effect is not for Pjj = (inverted 
hierarchy). The dotted curve represents the cases where the supernova matter effect is included 
but the Earth matter effect is not for Pjj = 0,1 (normal hierarchy) or Pjj = 1 (inverted hierarchy). 
The dot-dashed (dot-dot-dashed) curve represents the flux given by Eq. © (Eq. Q12|)) where the 
neutrinos do not oscillate at all (oscillate as in vacuum). 



